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Abstract. - We study the effect of the scattering of gravitational waves on planetary motions, 
say the motion of the Moon around the Earth. Though this effect has a negligible influence 
on dissipation, it dominates fluctuations and the associated decoherence mechanism, due to 
the very high effective temperature of the background of gravitational waves in our galactic 
environment. 



Decoherence plays an important role in the transition between microscopic and macro- 
scopic physics since it kills quantum coherences on a time scale which becomes extremely short 
for systems with a large degree of classicality . The details of this transition depend on 
the coupling mechanisms between the system under consideration and its environment. A lot 
of different models have been considered theoretically and decoherence has been experimen- 
tally observed in mesoscopic systems, such as a few microwave photons in a high-Q cavity S, 
for which the decoherence time is neither too short nor too long. 

For large macroscopic masses, say the Moon orbiting around the Earth, decoherence is 
so efficient that it precludes the observation of quantum coherences. It remains however 
important from a conceptual point of view to study the dominant decoherence mechanisms 
and to obtain a reliable estimation of the decoherence time scale. For motions in the solar 
system, decoherence is often attributed to the scattering of the electromagnetic fluctuations 
associated with solar radiation or cosmic microwave background. In this letter we show that 
the decoherence of planetary motions is not dominated by electromagnetic processes but rather 
by the scattering of stochastic gravitational waves present in our galactic environment. 

Gravitational fluctuations have already been proposed as a universal mechanism able to 
explain the transition from quantum to classical physics. If these fluctuations are characterized 
by length scales of the order of the Planck length ^,0, microscopic and macroscopic regions 

may be delineated by comparing the Planck length lp = £p and the Compton length 

lc = — associated with typical quantum phenomena for a mass m [js|^10[ . Remarkably, this 
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leads to a frontier determined by the Planck mass mp = y ^f, i-e. the mass scale built on 
the Planck constant h, the velocity of light c and the Newton constant G with a value 22 fig 
lying on the borderland between microscopic and macroscopic masses. However, this simple 
scaling argument is not by itself sufficient to reach any precise conclusion. 

In order to compare gravitational and electrodynamical contributions to decoherence, one 
has to discuss the corresponding fluctuation levels as well as their effects on the system, at 
the frequencies of interest for the latter. In this letter, we present a quantitative study of 
decoherence of planetary motions in the stochastic background of gravitational waves in our 
environment pl|,p^|. We give a precise estimate of the associated gravitational decoherence 
rate and show that it largely dominates the electromagnetic contribution. 

The interaction of macroscopic motions with gravitational fluctuations is well known from 
the theory of gravitational wave emission and gravitational wave detection [p~3| — p~5|] . At the 
limit of non relativistic velocities, the gravitational perturbation on a planetary system may 
be represented as a tidal acceleration acting on each mass 

x'{ (t) = -R m (t) Xj (t) R m (t) = -h^ (t) (1) 

The prime denotes a time derivative. The tidal tensor i?iojo is built up from components of 
the Ricmann curvature tensor with the index representing time components and indices 
i,j representing spatial components. -R;ojo has been written as the second order derivative of 
the metric perturbation evaluated at the center of mass of the planetary system in the 
transverse traceless (TT) gauge. The interaction is equivalently described by a Lagrangian 
perturbation coupling the quadrupole momentum of the system to the tidal tensor 

S'(t) = \h'( ] (t)Q i] (t) Q ii (t)=mL(t)x i (t)- 5 j-x k (t)x k (t) S j (2) 

The quadrupole Qij is reduced to its traceless part with a spatial Kronecker symbol. S' 
is the time derivative of the action integral as it is perturbed along a given trajectory, for 
example a circular orbit in the plane x%x 2 

Xl (t) = p cos (fit + 9) x 2 (t) = p sin (fit + 9) p 3 fl 2 = GM (3) 

For a system of two masses m a and m&, m denotes the reduced mass ™ a 7!£ lb , Xi the relative 
position and p the distance between the two masses. The last relation in (||) is the Kepler law 
which connects the radius p, the orbital frequency fl and the total mass M — m a + mj,. 

Our aim is to evaluate decoherence between two neighbouring motions on the same circular 
orbit. These motions correspond to slightly different values of 9 with the difference denoted 
by A9. The differential perturbation between the two motions is thus described by a quantity 
AS" written as in (gj) with Qy replaced by the difference AQy of the quadrupoles on the two 
motions. We write it as the product of a distance Ax and a force F 

1 x > (t) 

AS' (t) — F (t) Ax Ax = pA9 F (t) = -h'{x (t) mx\ (t) -^r- (4) 

2 J pfl 

Ax is the distance between the 2 motions, constant on a circular orbit, and F the component 
of the relative force projected along the mean motion. F may be expressed in terms of the 
circularly polarized metric perturbation h which fits the circular motion of the planetary 
system 



F(t) = ^= (h"e^ + h* "e-™*) h(t) = ^= [h 12 + (5) 
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The circular polarization h is normalized so that it corresponds to the same noise level as the 
two linear polarizations h\2 and h22 ~ hl1 in the case of an unpolarized background. Polariza- 
tions are defined with respect to the plane of the orbit and the two motions chosen in the 
vicinity of 9 = 0. 

We now come to a Fourier representation of the force F and metric h. Accordingly, force 
and metric fluctuations are described by noise spectra Cff and Chh 



F[uj} = J dt F (t) e luJt Cff [u] = J dt {F (t) ■ F (0)) e Mt 



(6) 



The dot specifies a symmetrical ordering when quantum fluctuations are considered Pqp7|| . 
The force (||) is driven by gravitational waves through a frequency transposition due to the 
evolution of the quadrupole momentum at frequencies ±20 

(u + 20) 2 h\uj + 20] + 20) 2 h* [uj - 20] 
F [u] = -mp ^-j= (7) 

We are mainly interested in the long term cumulative effect of fluctuations, that is in the 
momentum perturbation p t integrated over an interaction time t longer than the correlation 
time. The variance of p t is determined by the noise spectrum Cff [0] evaluated at zero 
frequency and expressed in terms of a momentum diffusion coefficient D 

Pt= [ dsF (s) (p 2 t ) = Cff [0] t = 2Dt (8) 
Jo 

The diffusion coefficient D gr due to gravitational waves is finally obtained as 

2D gT = Cff [0] = 4m 2 a 2 C hh [20] a = pQ 2 (9) 

a is the acceleration on the circular orbit and Chh [2fl] is the gravitational wave spectrum at 
frequency 2Q of the circular polarization h coupled to the system. Note that a gravitational 
background of galactic origin is certainly not isotropic whereas an extragalactic background 
may probably be considered as isotropic. For simplicity, forthcoming discussions are phrased 
in terms of an unpolarized and isotropic background but the general case is recovered by 
coming back to equation (||). 

The frequency of interest for the planetary motion of the Moon is |^ ~ 0.8 x 10~ 6 Hz. 
Information on stochastic gravitational waves around this frequency may be deduced from 
studies devoted to the detectability of gravitational background by interferometers [^8[|l9|| . 
We express the spectrum Chh as a gravitational wave energy /ceTgr or, equivalently, as a 
number n gr of gravitons per mode 



Chh [uj] = ^-k B T gT [u] = ^-g- Q + n gT [uj] ) hu> ( :i ) 



&b is the Boltzmann constant and T gr is an effective temperature of the gravitational back- 
ground. Figure 1 of El] and equation (3.1) of EJ| lead respectively to Chh [20] ~ 10 _34,5 Hz _1 
and Chh [20] ~ lO^'^Tlz -1 for |^ ~ 10 _6 Hz. Taking as a conservative estimate Chh — 
10 _34 Hz _1 , we obtain an effective temperature T gr w 10 41 K or, equivalently, a graviton num- 
ber per mode 10 57 . These numbers correspond to the high temperature limit 
fceTgr 3* Hence the gravitational noise is much larger than vacuum fluctuations which 
correspond to the term ^ in ( fio| ) and lead to ultimate fluctuations of geodesic distances of the 
order of Planck length |17| . 
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The estimations of Chh used here correspond to the confusion background of gravitational 
waves emitted by binary systems in our galaxy. They rely on the laws of physics and astro- 
physics as they are known in our local celestial environment. There also exist more speculative 
predictions for gravitational backgrounds emitted by cosmic processes [^| ^| . Depending on 
the parameters used in the models, these cosmic backgrounds may surpass the binary confu- 
sion background in the /iHz frequency range. Hence the latter can be considered as a minimum 
noise level in our gravitational environment. 

The diffusion coefficient may be written under the form of an Einstein fluctuation-dissipation 
relation, i.e. as the product of the effective temperature by a damping factor 

32Gma 2 . . 

D gI = mr gr fe B T gr r gr = 5c5 (11) 

The damping rate T gr is the inverse of the damping time associated with the emission of 
gravitational waves by the planetary system Jl3| — |l5|] . It does not depend on temperature and 
is so small, r gr m 10~ 34 s -1 for Moon, that it does not affect the classical motion on the orbit. 
Gravitational damping however has a noticeable effect in strongly bound binary systems such 
as millisecond pulsars [plj . 

At this point it is worth comparing the effects of gravitational and electromagnetic scat- 
tering. Modelling the moving mass as a sphere which perfectly scatters thermal photons at 
temperature T em , the damping rate r em associated with radiation pressure of these electro- 
magnetic fluctuations is evaluated as [E3| 



r..,.-^(*£=y (12) 

45m \ he J y ' 



cm 



The radius r has been supposed to be large compared to the wavelength of the photons. 
For the Moon, the temperature T om = 2.7K of the cosmic microwave background is already 
sufficient to produce a damping rate r cm «2x 10 _32 s _1 which is 200 times larger than r gr . 
Since r em varies as T e 4 m in agreement with Stcfan-Boltzmann law, the damping due to solar 
radiation is more than 10 10 times larger than r gr . In fact, the damping of the Moon, as it 
is revealed by the secular variation of lunar rotation p4j , is mainly due to the interaction 
between Earth and Moon tides and it corresponds to a damping rate more than 10 16 times 
larger than r gr . The role of gravitational scattering on damping of the Moon may thus be 
ignored but this is no longer the case for decoherence, as is shown in next paragraphs. 

In order to evaluate decoherence rates, we consider the effect of gravitational perturbation 
on the action difference ASt after an interaction time t 

AS t = f dsF (s) Aa; = p t Ax (AS 2 ) = (p 2 t ) Ax 2 (13) 
Jo 

The decoherence between the two neighbouring trajectories is measured by the mean value 
(exp lA ^* ) of the exponentiated dephasing. Since ASt is linearly driven by gravitational 
waves, it can be treated as a gaussian classical stochastic variable which leads to 

iASA ( (AS 2 t )\ , a 2 \ 

exp — — ) = exp I — -g- I = cxp (-A gr Aa; t) 

D RI 32Gm 2 a 2 , , , 

= "# = ^T k ^ (14) 
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Decoherence has been evaluated in the long term limit, where the force fluctuations may be 
approximated as a white noise characterized by a momentum diffusion coefficient D gr . This 
expression can also be derived by evaluating the Feynman- Vernon influence functionals p5| . 
often used in the study of decoherence p6[ , at the limit of high temperature. Note that 
the decoherence rate A gr is proportional to the square of the acceleration and, hence, would 
vanish if evaluated for an inertial motion. This has a simple interpretation through the 
Einstein formula (pT|). The damping rate r gr associated with the emission of gravitational 
waves indeed vanishes for inertial motion while the other factors entering the expression of 
-D gr do not depend on the specific motion. 

Using the Einstein fluctuation-dissipation relation ( jlT| ) and the expression (|l4|) of decoher- 
ence rates, the ratio between gravitational and electromagnetic contributions to decoherence 
can be rewritten 

A D V T 

A D V T 

ii em ^em x cm 1 em 

For the motion of the Moon, the gravitational decoherence rate A gr is found to be 10 38 larger 
than the rate A om corresponding to scattering of the cosmic photon background. It is still 
enormously larger than the effect corresponding to the scattering of solar photons. The same 
conclusion is reached for the comparison with the effect of tides. The latter effect determines 
the damping of the main motion of Moon but its contribution to decoherence is overshad- 
owed by gravitational scattering owe to the very high value of the effective gravitational 
temperature. As a consequence of this high temperature, the decoherence time is exceed- 
ingly short even for ultrasmall distances Ax. To fix ideas, the gravitational decoherence rate 
A gr « 10 75 s" 1 m~ 2 obtained for the motion of the Moon corresponds to a decoherence time in 
the 10/is range for a distance Ax between two trajectories as small as Planck length. 

The large value of gravitational temperature implies that decoherence should still be dom- 
inated by gravitational scattering for smaller size planetary systems. This can be discussed by 
writing the ratio between gravitational and electromagnetic decoherence rates as a product of 
dimcnsionlcss factors 



A gr _ 72 m 2 p 2 ( hfl \ 4 T gr 



A em 7r 3 mp r 2 \kBT cm J T ( 



(16) 



The factor ^ is clearly reminiscent of the simple scaling arguments presented in the In- 
troduction, illustrating the role of Planck mass as a reference on the borderland between 
microscopic and macroscopic masses. However the whole formula shows that these scaling 

2 

arguments are not sufficient for obtaining reliable quantitative predictions. The factor ^ is 
a geometrical factor depending on the radius p of the orbit and the radius r of the orbiter. 
Then the ratio ( |l6| ) also depends on the inverse fourth power of the photon number ^J™ 1 
per mode at the orbital frequency and on the already discussed ratio between temperatures 
of graviton background and photon background. Clearly the last three dimensionless factors 
appearing in ( |l6| ) have no relation with the scaling arguments of the Introduction but have to 
be known for a quantitative comparison of gravitational and electromagnetic contributions. 
If we consider as an example a man made gravitationally bound planetary system consisting 
of two spheres having an ordinary metallic density, we obtain comparable decoherence rates 
for gravitational scattering and electromagnetic scattering of cosmic microwave background 
in the case of a compact geometry with masses of the order of 10 3 kg. 

We have shown in this letter that the scattering of gravitational waves in our celestial 
environment is a dominant cause of decoherence for planetary motions such as the motion 
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of the Moon around the Earth. Due to the very large effective temperature of the gravita- 
tional wave background, this mechanism leads to a decoherence by far more efficient than 
the other fluctuation mechanisms though its contribution to damping of the mean motion 
can be neglected. As far as the theory of measurement is concerned, this implies that the 
ultimate fluctuations of the motion of Moon are determined by the same classical gravitation 
theory which also explains its mean motion. Precisely these fluctuations are determined by 
the classical gravitational waves present in our local celestial environment. 

* * * 
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